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A detailed study of the anharmonicity-induced resonances caused by the coupling of center-of- 
mass and relative motion is presented for a system of two ultracold atoms in single-well potentials. 

As has been confirmed experimentally, these inelastic confinement-induced resonances are of interest, 
since they can lead to coherent molecule formation, losses, and heating in ultracold atomic gases. 

A perturbative model is introduced to describe the resonance positions and the coupling strengths. 

The validity of the model and the behavior of the resonances for different conhnement geometries are 
analyzed in comparison with exact numerical ab initio calculations. While such resonances have so 
far only been detected for large positive values of the s-wave scattering length, it is found that they 
are present also for negative s-wave scattering lengths, i. e. for attractive interactions. The possibility 
to coherently tune the resonances by a variation of the external confinement geometry might pave 
the way for coherent molecule association where magnetic Feshbach resonances are inaccessible. 


I. INTRODUCTION 

Theoretical treatments of strongly-correlated ultracold 
atoms in single-well potentials routinely adopt the har¬ 
monic approximation to describe the trapping potential. 
This is, of course, an idealization, because every realistic 
trapping potential is finite and thus anharmonic. A ma¬ 
jor benefit of a harmonic confinement is that for identi¬ 
cal particles relative (rel.) and center-of-mass (c.m.) mo¬ 
tion decouple. Moreover, especially for deep potentials 
the harmonic confinement resulting from a second-order 
Taylor expansion of the potential in its minimum might 
be a good approximation, especially since the center-of- 
mass to relative motion (c.m.-rel.) coupling introduced 
by the anharmonicity of the trapping potential is ener¬ 
getically negligible compared to the energy scale of the 
interatomic interaction. As a consequence, the theoret¬ 
ically predicted binding energy of two ultracold atoms 
confined in a harmonic trap [1] has been experimentally 
confirmed [5]. Moreover, the measurement [3] and cal¬ 
culation [3] of the influence of the anharmonicity on the 
energy of states in deep optical lattices has revealed that 
the deviation to the harmonic approximation for the low¬ 
est band is negligible in most cases. 

On the first glance, it was therefore surprising that 
c.m.-rel. coupling can have a significant impact on an ul¬ 
tracold atomic quantum gas. In mm it was revealed that 
the particle loss and heating observed in [7] was caused 
by inelastic confinement-induced resonances (CIR), i. e. 
resonances due to the c.m.-rel. coupling of a c.m. ex¬ 
cited bound state with an unbound atom pair - even in 
a deep optical lattice where the harmonic approximation 
has proven accurate for the rel. motion energies. The ex¬ 
planation of the losses due to inelastic CIR was necessary 
in order to resolve contradicting results of the experiment 
[7] and the theory of elastic CIR [H E] that originally 
was used as an explanation for the observed atom losses. 
Additionally, alternative explanations of the experimen¬ 
tal findings were proposed based on the assumption of a 
harmonic confinement. One is based on multichannel ef¬ 


fects m, others [3 HI] on a Feshbach-type mechanism. 
Finally, in |12j it was demonstrated by simultaneously 
observing both the elastic and the inelastic CIRs and by 
excluding other mechanisms like many-body effects due 
to the design of the experiment that a coherent molecule 
formation is triggered uniquely at a c.m.-rel. coupling 
resonance which confirms the explanation of the losses in 
[7| in terms of the inelastic CIR. This demonstrates the 
importance of the anharmonicity-induced c.m.-rel. cou¬ 
pling, especially in view of its universality. 

In fact, very recently it was demonstrated that inelas¬ 
tic CIR are also present in ultracold dipolar systems m, 
and even for Coulomb-interacting systems such as exci- 
tons in quantum-dot systems M- For the latter, the 
resonances were proposed to serve as a novel kind of 
controlled single-photon source. At the inelastic CIR a 
variation of the exciton confinement leads to a redistri¬ 
bution of the charge density with subsequently increased 
annihilation probability of the electron-hole pair. Since 
this process can be steered in situ by a variation of the 
external confinement, single photons can be emitted on 
demand. 

Resonances in atomic gases due to c.m.-rel. coupling 
have been mentioned in literature before. The first work 
explicitly discussing a possible molecule formation due to 
anharmonicity-induced c.m.-rel. coupling is m- In the 
model in [T5| a deep isotropic optical lattice is consid¬ 
ered. The evaluation of the c.m.-rel. coupling matrix el¬ 
ements is performed with wavefunctions in the harmonic 
approximation. A direct loss process is considered where 
two unbound atoms couple to a molecule in the contin¬ 
uum, i. e. the c.m.-part of the bound state is assumed to 
be a highly excited c.m. wavefunction which is approx¬ 
imated by a spherical wave. The work concludes that 
the dimer-production rate is too small to be signihcant 
in an optical-lattice experiment. The observed losses in 
[7| proof different, however, i. e. losses at inelastic CIR in 
a deep optical lattice are observed. The reason for the 
small dimer-production rate in |I5j is that the coupling 
to highly excited bound states is very small as will be 
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shown in this work. On the other hand, as demonstrated 
in [ana, losses at inelastic CIR occur dominantly due to 
the coupling of a rel.-motion bound state with low-order 
c.m. excitation and subsequent three-body collisions, in 
contrast to a direct coupling to a rel.-motion bound state 
in the c.m. continuum as assumed in [Ill- 

Even in a harmonic confinement a coupling of the c.m. 
and rel. motion is present in the case of heteronuclear 
atoms or identical atoms in different electronic states 
[MU]- Moreover, the occurrence of Feshbach-type res¬ 
onances due to the c.m.-rel. coupling was discussed in 
m, their behavior in a superlattice was characterized 
in m- In mixed dimensions, the experiment performed 
in m also detected inelastic loss resonances for a varia¬ 
tion of the scattering length. The behavior of the elastic 
and inelastic CIR within a quasi-ID lattice model was 
considered in |22) . 


In the present work, the two-body model that was in¬ 
troduced in [S] for the description of the positions of 
the resonances in quasi-lD and quasi-2D confinement ob¬ 
served in [7] is generalized to describe the position and 
coupling strength of in principle all c.m.-rel. coupling res¬ 
onances in single-well potentials of arbitrary symmetry. 
The validity of the model is discussed in comparison to 
ab initio calculations. Resonance positions and coupling 
strengths are investigated for a change in the confine¬ 
ment geometry, i. e. in the transition between an almost 
isotropic to a cigar-shaped (quasi-ID) confinement and in 
the transition between an almost isotropic to a pancake¬ 
shaped (quasi-2D) confinement. It is demonstrated that 
higher-order resonances occur also for negative values of 
the s-wave scattering length. Consequently, by analyzing 
the wavefunction of the system it will be demonstrated 
that molecule formation also occurs in the strongly at¬ 
tractive interaction regime. 


The paper is organized as follows: In Section [IT] the 
two-body Hamiltonian is introduced and the basic idea 
of inelastic CIR is briefly recapitulated. Then a model 
is formulated to predict th e p osition (Section |III[ ) and 
coupling strength (Section IVI of c.m.-rel. resonances. 
To validate the model, the behavior of the resonances is 
discussed under a varying geometry of the confining po¬ 
tential in Sections El and EB Section IVIII deals with the 
optimization of c.m.-rel. coupling and the limitations of 
the model. It is then concluded in Section ElIIl that res¬ 
onances ignored in previous considerations (e. g. in [5]) 
should also lead to a molecule formation in the strongly 
attractive interaction regime of ultracold atoms. The ar¬ 
ticle closes with a conclusion in Section [ 





Figure 1. (Color online) Eigenenergy spectrum of two ultra¬ 
cold atoms in an isotropic harmonic trap interacting via a 5 
pseudopotential PUS] for a varying s-wave scattering length 
a (or varying confinement length dho which would, however, 
also change oj continuously). The upper panel shows the en¬ 
ergy of the rel. motion Hamiltonian H^ei, the lower one the 
combined energy of the rel. and the c.m. motion Hamiltoni¬ 
ans, i. e. idrei + Hem- Introducing a coupling between the rel. 
and c.m. motion makes the crossings avoided as illustrated 
in the lower part as solid lines, while the black dashed lines 
indicate the diabatic curves. Passing through the crossing 
adiabatically (on the solid line) allows for a transformation of 
the bound state with c.m. excitation into a trap state in the 
c.m. ground state. 


II. THE HAMILTONIAN 

In c.m. and rel. coordinates, r = ri — r 2 and R = 
i(ri -|-r 2 ), respectively, the Hamiltonian of two identical 
particles in an external trapping potential can be written 


R(r, R) = Rrei(r) + Rcm(R) + W(r, R) (1) 
Rrei(r) = T,ei(r) + Vei(r) + 17i„t(r) (2) 

Rcm(R) = Te„,(R) + Vcm(R) ■ (3) 
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where Trei and Tcm denote the kinetic-energy operators 
of the rel. and c.m. motion, respectively, and Vrei and Vcm 
are the separable parts of the potential energy. Thus, W 
contains only the non-separable terms of the potential 
energy. Uint{r) is the interparticle interaction. 

In Fig. 2 the eigenenergy spectrum of the Hamiltonian 
in Eq. (IT I for two ultracold atoms interacting via the S 
pseudopotential 


Uint{r) 


Airfi^a 

m 


6{r) 


dr 


(4) 


a being the s-wave scattering length and m the atomic 
mass, is shown. The atoms are confined in an isotropic, 
i. e. spherically symmetric, harmonic trapping potential. 
The energies are plotted for a varying inverse s-wave scat¬ 
tering length dho/a, being the harmonic os¬ 

cillator length. The spectrum of iJrei (Eq- (§) contains 
a bound state \ipb) bending down to negative infinity for 
a —>■ O’*" and trap states, the energetically lowest one de¬ 
noted as IV’i)- In case of a harmonic trapping potential, 
the coupling W between the c.m. and the rel. motion van¬ 
ishes. Hence, in order to obtain the spectrum of the full 
Hamiltonian (Eq. 0 ) the c.m. energies are added to each 
state of the rel. spectrum, resulting in the middle plot of 
Fig-B Each state of the rel. spectrum appears now with 
an infinite series of c.m. excitations. Crossings appear 
between c.m. excited bound states, e. g., { 2 , 0 , 0 })^ and 

trap states, e.g., the lowest trap state |^/’i‘i>(o,o,o))- In¬ 
troducing a coupling VF(r,R) 7 ^ 0 between c.m. and rel. 
motion turns the crossings to avoided ones (besides mod¬ 
ifying the energies of the states also globally) and enables 
an adiabatic transition of a c.m.-excited molecular state 
and a state of an unbound atom pair in the c.m. ground 
state a indicated in the lower plot of Fig. B 

Therefore, the c.m.-rel. coupling introduces a 
Feshbach-type resonance at the crossing position. The 
occupation of the bound state at the resonance is only 
possible because the excess binding energy can be 
transferred into c.m. excitation energy due to the anhar- 
monicity of the confining potential. This redistribution 
of binding energy to kinetic energy marks an inelastic 
process and thus these c.m.-rel. coupling resonances were 
denoted as inelastic CIR. In m it is demonstrated how 
a coherent molecule formation can be realized at the 
c.m.-rel. coupling resonance. Such a molecule formation 
can trigger particle loss and heating in a many-body 
system in a two-step process. At the resonance two 
atoms coherently couple to the c.m.-excited molecular 
state. Then, this molecule collides either with another 
molecule or an unbound atom leading to a deexcitation 
of the molecule into a deeply bound state and subsequent 
loss of the involved particles from the trap. 

The full two-body spectrum Fig. B shows a plethora of 
crossings. A full six-dimensional treatment of the two- 
body problem involving the c.m.-rel. coupling is possible 
in mil but numerically quite demanding. Since c.m.-rel. 
coupling resonances can have a substantial influence on 


the stability of an ultracold atomic gas, the knowledge 
of the position and coupling strength of the resonances 
is of great interest. The understanding and assignment 
of these resonances is in particularly important for an 
unbiased identification of other few-body resonances of, 
e. g., Efimov type. Therefore, simplified models that al¬ 
low for the estimate of resonance positions and widths 
are desirable. 


III. C.M.-REL. COUPLING MODEL - 
RESONANCE POSITIONS 


A model is introduced to predict the position and cou¬ 
pling strength of the inelastic CIR. The theory general¬ 
izes the model described in [S] that is only applicable to 
strongly anisotropic (quasi-ID or quasi-2D) confinement 
and does not consider coupling strengths explicitly ^5] . 

Of course, for c.m.-rel. coupling resonances to be 
present a c.m.-rel. coupling must be introduced. There¬ 
fore, the harmonic approximation has to be abandoned. 
Optical-lattice potentials [55] are widely used in ultra¬ 
cold experiments and offer a great degree of flexibility 
and control. In a deep optical lattice, i. e. ^ 1 where 

t2 J.2 

Er = is the recoil energy and V is the lattice depth, 
tunneling between neighboring wells is suppressed and 
the potential can be regarded as a stack of single-well 
potentials. In this work, a sextic potential is used re¬ 
sulting from a Taylor expansion of a sin^ optical-lattice 
potential up to the sixth degree. A separation of the 
expansion in rel., c.m., and coupling terms 


ITei(r)= ^ U, 

j^x,y,z 

Icm(R-) = Vj 

W{r,R)= 


1/^2 2 _ J_l4 4 , J_;U6 6 

2 ^ ^ 24 ^ ^ ^ 720 ^ ^ 


(5) 


+ ( 6 ) 

(7) 


respectively, shows that the quartic terms all have a nega¬ 
tive sign which makes the expansion to the sextic degree 
necessary. Otherwise, an unphysical continuum would 
occur in the spectrum reaching in energy towards neg¬ 
ative infinity. In Eqs. (B)-(B) Vj is the lattice depth in 
direction j G {x, y, z}, kj = and Xj is the wavelength 
of the laser in direction j. Introducing the harmonic os¬ 
cillator frequencies oJj = the potential terms can 

be written in a more canonical form. 

It has been demonstrated mum that sextic poten¬ 
tials are well suited to describe anharmonicity-induced 
c.m.-rel. coupling in single-well potentials. The large 
flexibility in the potential parameters makes the sextic 
potential, moreover, suitable to accurately describe the 
c.m.-rel. coupling in a variety of potentials, e. g. Gaussian 
beam potentials m- 
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It is first assumed that the anharmonicity only has a 
small (global) influence on the eigenenergies of the states. 
Hence, the harmonic approximation is used and the po¬ 
sition of a c.m.-rel. coupling resonance is determined by 
the position of the crossing of the c.m. excited bound 
state \ipb^n), where n = {nx,ny,nz) is the quantum 
number of the c.m. excitation, that separates spatially 
for an optical-lattice potential, and a trap state 
as illustrated in Fig. 

The energy SJ'®* of the bound state in an harmonic 
confinement of arbitrary anisotropy in dependence of the 
s-wave scattering length a is given implicitly by [27] 



a 



\/(l - e-‘) (1 - ‘) (1 - *) 



dt 


( 8 ) 

where e = (Ff),'®' - Eo)/{Hujy), Eq = + cvy + w^), 

Vx = and r]z = ujziujy. The difference to other 

equations for the bound state, like e. g. in [1] (valid for 
a 3D isotropic confinement), in (valid for a 3D har¬ 
monic confinement of single anisotropy), or in [28j (valid 
for only transversally trapped atoms), is that it is valid 
for an arbitrarily anisotropic 3D confinement. 

A general expression for the eigenenergies of the trap 
states in an arbitrarily anisotropic confinement, i. e. 
states above Eth = ^{oJx+^y+^z), is not known yet. For 
ultracold temperatures the occupation of excited states 
is suppressed. Hence, in the following only crossings with 
the first trap state |^i) are considered. Assuming with¬ 
out loss of generality that mm{ujx,ujy,ujz) = uJz (unless 
stated differently), the eigenenergy A)®' of |V'i) lies in 
the interval [Eth, Eth + 2hujz)- It can be shown that in 
the case of an isotropic harmonic confinement the cross¬ 
ing between an excited bound state with a single lowest- 
order c.m. excitation (n^ = 2,nj^i = 0) with the first 
trap state occurs at 


i?)®^ — Eth + JEJz (9) 

which is thus chosen for the model as an approximation 
of the energy of the first trap state. For crossings with 
higher trap states, the model can be extended by the 
proper choice of the energy of that trap state. 

For a spatially decoupled confinement, like ex¬ 
pansions of an optical-lattice potential, the eigen¬ 
states of the c.m. Hamiltonian factorize as $„(R) = 
</'nJ-A)())„^(y)())„jZ) with eigenenergies E^ = 
^j=x y z ^3 + 5 )- When combining rel. and c.m. mo¬ 

tions the energies of the bound states become El^^{a) -I- 
A®™ while the energy of the lowest trap state is given by 
A)®' -|- Q). Crossings between a c.m. excited bound 
state and the lowest trap state are determined by 


( 10 ) 


where 


A _ jccni iccm 


( 11 ) 


is the c.m. excitation. The corresponding s-wave scatter¬ 
ing length at the crossing is obtained from Eq. Q. 

So far all energies were treated within the harmonic 
approximation. It will be demonstrated that this purely 
harmonic model gives in some cases already good quan¬ 
titative results. However, for the c.m. excitations A„ 
higher c.m. states are involved. Additionally, for small 
A„ the bound state crosses the trap state for small pos¬ 
itive dho/cL or even negative dho/o- where the two states 
cross with comparable slopes, see Fig. [l] For such cross¬ 
ings the position is very sensitive to the energies of the 
involved states. Hence, the energy of the first trap state 
as well as the c.m. excitation A must be corrected. 


Treating the anharmonic terms of the c.m. potential in 

^ ^ i?® (here written 


Eq. n, E3=x,y.z-T,^R^ + J20 V: 


•j=x,y,z 24 Vj ""3 720 3 

3 

in dimensionless units of energies in Hujj and lengths in 
with M = 2m) within first-order perturbation 
theory results in 


A, 


= Y. 


3=^,V,z 


hujj 


(n/ +nj) 


fe2 2 

576y2 + 3n/ +4nj) 


( 12 ) 


for the c.m. excitation and 


EJ®' = Fw, + 


E 

j^x,y,z 


1 ^ hujj hiOj 


for the energy of the first trap state (see Section XB 
the Appendix for details). 

It will be demonstrated that in the case of resonances 
for negative values of the s-wave scattering length and 
strongly anisotropic confinement, the effective ID c.m. 
problem needs even to be solved exactly and thus nu¬ 
merically to obtain accurate results. For the numerical 
evaluation of the stationary ID Schrodinger equation, the 
approach described in was used. 

Therefore, three models for the resonance position were 
introduced that differ by the treatment of the c.m. exci¬ 
tation A and by the energy of the first trap state E)®b In 
model A, A and_K®' are given in the harmonic approxi¬ 
mation by Eqs. ([TT| and Q, in model B by Eqs. ( 12 ) and 


(131 within a perturbative correction, respectively, and in 


model C the c.m. energies are calculated numerically ex¬ 
act. For given A and E®®* within one model, the inelastic 
CIR position for a c.m. excitation is then determined by 
Eq. (pj). 


^rel ^ ^rel _ ^ 
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IV. C.M.-REL. COUPLING MODEL - 
COUPLING STRENGTHS 

After having introduced a straightforward procedure 
to evaluate the resonance positions using a model (with 
three different versions A, B, and C) for the resonance po¬ 
sitions, the coupling strengths are considered. They are 
of particular interest for experiments, since they deter¬ 
mine the width and thus experimental visibility. Further¬ 
more, they are, e. g., necessary for a Landau-Zener treat¬ 
ment of the dynamics at the resonances that allows for 
an estimate whether diabatic or adiabatic transitions be¬ 
tween the involved states occur at the avoided level cross¬ 
ings. The matrix element defining the coupling strength 
Wn between a bound state <i>n(R)) with c.m. ex¬ 

citation An and the lowest trap state |i/'i(r) $(q o,o)(R')) 
is 


1T„ = (V^('’)(r)<i>„(R)|IF(r,R)|V^i(r)$(o,o.o)(R)) • 

(14) 


For the model the wavefunctions of a harmonic confine¬ 
ment are adopted. Hence, the c.m. wavefunction is the 
product of ID harmonic oscillator wavefunctions (here 
written in dimensionless units of energies in Hujj and 


lengths in 



$„(i?j) = TT ■! 


‘"nl 


e-!"^ RniRj) 


(15) 


function, and + y^. The wavefunctions are writ¬ 

ten in dimensionless units of energies in hujz and lengths 

in dz = Moreover, the previously introduced def¬ 
initions ^ and e = — | {lOx + ojy + uJz) remain 

valid. Within the model (see above), the energy at the 
resonance is e = Hujz (assuming that the elongation of 
the trap is along the z direction). It is important that 
for the wavefunction the energy is not corrected by the 
anharmonic terms, because this might result in a differ¬ 
ent energy branch of the spectrum. With e = hujz at the 
resonance, Eq. 0 can be further simplified (in physical 
units) to 


tpip^z) 


V 2 r(-|) 

47r2d5_-\/dz 


\z\ exp 




(18) 


with d± = y and using that for integer n the rela¬ 
tion D„(a;) = (^) 

In quasi 2D, the wavefunction can be simplified to |23j 

V. (A *) = (- i) ^, 1 , ) 

(19) 


where U denotes the confluent hypergeometric function. 
Again, dimensionless units of energy in HlUz and lengths 
in dz are used. 

The rel. motion bound-state wavefunction [23] 


where Hn(i?j) denote the Hermite polynomials. The c.m. 
integral in Eq. (141 reduces to a sum of ID integrals that 


can be calculated even analytically. 

The 3D integral over the relative-motion coordinate is 
more laborious. While an expression for the trap wave- 
functions for an arbitrarily anisotropic harmonic con¬ 
finement is so far (to the authors’ knowledge) not yet 
known, a general solution for the trap state in a har¬ 
monic potential with a single (but arbitrary) anisotropy, 
e. g. uJx = Wy =: a;_L ^ Wz is given in [23]. However, the 
numerical evaluation of Eq. (14) with the most general 


version of the relative motion wavefunction given in 


2^3/22.72 ^ 

OO 




m—0 


/ 2mrj — e 

V 2 


^e-2mr, , 

(16) 


has turned out prohibitively demanding. In the regime of 
a strongly elongated (quasi-lD) potential the expression 
greatly simplifies [53| to 

'■(-I) D. (v^lzl) . (17) 

In these equations denotes the parabolic cylinder 
function, the Laguerre polynomials, F the gamma 


i’bir) = 


74 


dj_^{2Trf^^ 


dt ■ 


exp( 


2tdz 


^coth 


(&J) 


4 sinh(|^) 


( 20 ) 


which is written here in in physical units, is valid for an 
isotropic confinement as well as for a strongly anisotropic 
trap geometry. Note, if the energy of the bound state 
crosses the energy of the trap state for positive values 
of the s-wave scattering length that are sufficiently small 
(a 7 '^-l)) its energy at the crossing is sufficiently small 
such that the bound state wavefunction can be described 
in good approximation by its trap-free counterpart 


1plree{r) = 


y/2jra 


( 21 ) 


However, for crossings at negative s-wave scattering 
lengths, this approximation certainly fails, because in the 
free case the bound state only exists for a > 0. Hence, in 
the following Eq. (201 is used for the bound-state wave- 
function. 

The wavefunctions of a harmonic confinement of sim¬ 
ple anisotropy are used for the model. Therefore, the 
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rel. motion integral is reduced to a two-dimensional one 
where the coupling term W^(r, R) is averaged over the 
transversal direction, uj± = ^{ujx +^y)- The matrix ele¬ 
ment becomes 

poo poo 

Wn = 27r / dp p dzil^b{p,z)ilJi{p,z)W{p,z) (22) 

0 J —oo 

with 


/ OO 

dRkWkip, 

K—x,y 

/ oo 

dZW,{z,Z)^n^<^>n, (23) 

-OO 

where VF(r,R) = ^j) n = 

( 77 , 3 ^, Tly^ ’^z)' 

Hence, different to the resonance position where three 
models were introduced that differed in the treatment of 
the c.m. energies, two models for the coupling strength 
were introduced that differ in the way the rel.-motion 
wavefunction is treated. In quasi ID, Eq. (18) is used 
for ijji within model 1, in quasi 2D Eq. (19) is adopted 
for ipi within model 2. The models solve the coupling 
matrix element of Eq. (|14| for the exact coupling term of 


the sextic potential of Eq. ([^ with wavefunctions in the 
harmonic approximation. 


Wx/wy = 1.1, Vy = 35.9 Er- To obtain an elongation in 
the longitudinal z direction, the potential depth Vz is de¬ 
creased. Hence, an almost spherical potential is deformed 
into an elongated, cigar-shaped one. 

The coupling term W in Eq. ([^ is totally symmetric 
and hence only states of equal symmetry couple. Since 
in the following only crossings with the first trap state 
without a c.m. excitation are considered, only even exci¬ 
tations, i.e. states with even n^, can couple. Moreover, 
in the ab initio calculation it suffices to consider the Ag 
symmetry because the states involved in the considered 
inelastic CIR, the first trap state, the last bound state 
and the even c.m. excitations possess Ag symmetry, see 
[24] for details. 


The two lowest-order transversal resonances with 
(2, 0, 0) and (0, 2, 0) c.m. excitation and the first-order 
longitudinal resonance with (0,0,4) c.m. excitation are 
investigated in the following. These resonances are se¬ 
lected because they are the most pronounced ones. As 
a rule of thumb, the coupling between a c.m. excited 
bound state and the ground trap state decreases with 
the order of the c.m. excitation, i. e. the lowest reso¬ 
nances Ui = 2,nj^i = 0 show the strongest coupling. 
A simple argument for this rule is that the stronger os¬ 
cillations in higher excited c.m. bound states decrease 
the overlap to the trap state and hence also the coupling 
matrix element Eq. (14). Numerically, this rule is veri- 


V. 3D-1D TRANSITION 


n 

(0,2,0) 

(0,2,2) 

(0,4,0) 

(0,6,0) 

OO 

0 

Wn[W~^hu}y] 

3.19 

1.46 

0.795 

0.50 

0.31 


In order to characterize the behavior of the reso¬ 
nances and to validate the models, full six-dimensional 
ab initio calculations are performed. The numerical al¬ 
gorithm to solve the stationary Schrodinger equation 
for the Hamiltonian in Eq. ([^ is described in |2^ . 
Eor the efficient computational treatment of ultracold 
atoms in optical lattices the basis functions are sym¬ 
metry adapted to the eight irreducible representations 
(Ag,Rig,i? 2 g,R 3 g, Au,i?iu,i? 2 u,R 3 u) of the orthorhom¬ 
bic point group D 2 \i. This leads to a corresponding block 
structure of the Hamiltonian matrix. As a realistic ex¬ 
ample for an interatomic interaction potential t/int the 
(numerically given) Born-Oppenheimer potential curve 
of two ^Li atoms in the electronic state a are used in 
this work, see m for details on how this potential was 
obtained. The variation of the s-wave scattering length 
is achieved by slight modifications of the inner-wall of 
the potential [30] which effectively changes the asymp¬ 
totic behavior of the radial wavefunction and hence the 
s-wave scattering length. Resonance positions and cou¬ 
pling strengths are extracted from the ab initio data by 
fitting a two-channel model to the corresponding avoided 
energy crossing (see Section X A in the Appendix). 


In the following, the behavior of the coupling reso¬ 
nances is investigated in the transition from a 3D to 
a quasi-lD confinement. Two ^Li atoms are consid¬ 
ered in a sextic potential (Eqs. ©-0) with A=1000nm, 


Table I. Coupling strengths Wn for uJxjuJy = 1.1, uJyloJz = 1, 
Vy = 35.9 Er of the first trap state with c.m. excited bound 
states for different c.m. excitations n obtained with the model 

1 . 


fied by calculating the coupling matrix element Eq. (14) 
for different c.m. excitations, see Table|l| confirming that 
Wn decreases with an increasing order of the c.m. excita¬ 
tion. It can also be seen in the energy spectrum in Fig. 
where high-order resonances do not show visible avoided 
crossings |3T|. The fact that the coupling decreases with 
the c.m. excitation of the bound state delivers an expla¬ 
nation why in m it is (correctly) concluded that the 
dimer-production rate at an c.m.-rel. coupling resonance 
involving a very highly c.m.-excited bound state is neg¬ 
ligible. 

While lowest-order resonances show the strongest cou¬ 
pling, in the following the lowest-order longitudinal reso¬ 
nance with c.m. excitation (0,0, 2) is not considered be¬ 
cause its position fades away to ^ <C — 1 with decreasing 
ujz where the bound state is very shallow and has lost its 
characteristic small interatomic distance. In the full spec¬ 
trum the crossing can therefore not be easily resolved any 
more. This can be seen in Fig.[^in the upper most panel, 
where the energy of the ( 0 . 0 , 2 )) state asymptotically 
approaches |'i/'i‘^’(o,o,o)) without a pronounced crossing. 
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Figure 2. (Color online) Ab initio energy spectrum of a 
^Li-^Li system confined in a sextic trapping potential with 
A=1000nm, uj^ltOy = 1.1, = 10, Vy = 35.9 i?r. The la¬ 

bels of the energy branches denote the corresponding diabatic 
states. 

In Fig. 1^ the three considered resonance positions of 
the models with different degrees of corrections together 
with ab initio results are shown. Before discussing the 
validity of the models, the behavior of the resonances is 
analyzed. For small anisotropies ooyjujz < 2 the (0,0,4) 
resonance lies at larger values of dy/a (the green curve 
lies above the red and blue curves in Fig. simply be¬ 
cause the c.m. excitation A(q q 4 ) is larger than A(q 2 , 0 ) 
and A( 2 _o,o)- Therefore, the c.m. excited bound state 
l'0ft‘J’(o,o,4)) crosses at larger values of dy/a than the ones 
with a single excitation. A decrease of I 4 oc decreases 
the spacings of the states that have a c.m. excitation in 
the z direction. Hence, for decreasing the (0,0,4) 
resonance crosses constantly at smaller values of dy/a 
which explains the monotonic decrease of the ab initio 
results in Fig. (green squares). On the other hand, 
the transversal c.m. excitations (A(o_ 2 ,o) and A( 2 .o,o)) re¬ 
main constant for a variation of Yet, a decrease in 
also decreases the energy of the first trap state. There¬ 
fore, the transversally excited c.m. states |'0b4>(2,o.o)) and 
l^b‘J’(o,2,o)) cross at larger values of dy/a with decreasing 
ujz converging to a finite value as —>■ 0. 

Next, the validity of the models is considered. For 
the transversal (2,0,0) and (0, 2,0) resonances the per¬ 
turbative corrections Eqs. ( |l^ and ( |l^ agree with the 
numerically exact corrections (the dashed and solid lines 
are indistinguishable). The resonance positions resulting 
from the harmonic approximation (dashed-dotted lines, 
model A) show an almost constant offset towards larger 
values compared to model B were the anharmonicity in 



Figure 3. (Color online) Resonance positions for different c.m. 
excitations n = (nx,ny,nz) (see legend) obtained by full ab 
initio calculations [dots (2,0,0), stars (0,2,0), squares (0,0,4)] 
and the model (lines) if the tightness in 2 direction is re¬ 
duced. The dashed-dotted lines correspond to the harmonic 
approximation (model A). The dashed lines are obtained by 
a perturbative correction of the energies (model B). For the 
solid lines, the correction is calculated in a numerically exact 
way (model C). 


the C.m. motion has been taken into account. This small 
offset is due to the missing negative quartic term that 
is present for the sextic potential. Certainly, the models 
give very good quantitative agreement with the ab initio 
calculations. For strong anisotropies the results of the 
models B and C (dashed and solid lines) are in perfect 
quantitative agreement, e. g., at ujy/wz = 10 the models 
give dy/a = 1.397 and the ab initio method results in 
dy/a = 1.396 for the (2,0,0) resonance. This excellent 
applicability of the model for the resonance positions in 
the quasi-lD regime lead to the quantitative agreement 
of the model compared to the loss resonances measured 
in [7] as shown in [^. 

The results for the longitudinal (0,0,4) resonances 
are more sensitive. First, with a decreasing potential 
depth the anharmonicity is important already for low 
lying states. Second, for the (0,0,4) resonance higher 
c.m. excitations are involved which enhances the influ¬ 
ence of the anharmonicity. Moreover, for a decreasing 
resonance position the bound state crosses the trap state 
with an increasingly comparable slope which makes the 
position more sensitive to energy variations. Therefore, 
model A including uncorrected, harmonic c.m. excitation 
is inaccurate over the entire range of anisotropies. The 
perturbatively corrected model B is acceptable for mild 
anisotropies (uiy/iOz ^ 10) but has a wrong behavior for 
ojy/ujz /Z 10. Finally, model C which corrects the c.m. 
excitations and the energy of the trap state exactly nu¬ 
merically is quantitatively accurate over the entire range 
of the scattering length, even in the limit uiz —>■ 0 (green 
solid line). 

The quantitative accuracy for very large values of the 






















anisotropy, ujy/uiz > 10, could not be confirmed by the 
ab initio method [ 21 ] • The basis set of the method con¬ 
sists of spherical harmonics which are not well suited 
for resolving extremely anisotropic wavefunctions unless 
high angular-momentum quantum numbers are employed 
which leads to an inconvenient scaling of the numerical 
effort. 



Figure 4. (Color online) Coupling strength for different c.m. 
excitations n = (n^;,n^,n^) (see legend) obtained by full ab 
initio calculations [dots (2,0,0), stars (0,2,0), squares (0,0,4)] 
and the model 1 (solid lines). 


In Fig.j^the coupling strengths corresponding to Fig.[^ 
are shown. Again, the overall behavior is discussed prior 
to the validity of the model 1. The coupling strength 
of the transversal resonances with (2,0,0) and (0,2,0) 
c.m. excitations (red circles and blue stars, respectively) 
decreases monotonically with increasing anisotropy but 
approaches a constant value for Wz —> oo. Otherwise the 
observation of the (2,0,0) and (0,2,0) resonances in [7] 
where the anisotropy of the external confinement is « 10^ 
would not have been possible. A simple argument for the 
monotonic decrease is that VF(r,R) oc 14. 

For the resonance with longitudinal (0,4, 0) c.m. exci¬ 
tation, a non-monotonic behavior is visible. In the limit 
of Wz —>■ 0 which corresponds to a zero potential in the 
z direction, the coupling of the resonances with a c.m. 
excitation in the longitudinal direction vanishes. This is 
intuitive, since without a confinement potential there ex¬ 
ists no confinement-induced c.m.-rel. coupling. On the 
other hand, a decrease in the potential depth 14 leads to 
an enhancement of the anharmonicity-induced coupling, 
since the potential becomes more anharmonic (this will 
be discussed in more detail in Section [VlT ). The result of 
these counter-acting effects is the non-monotonic curve 
with the local maximum for the (0,0,4) resonance and a 
vanishing coupling for Wz —f 0. 

Next, the validity of model 1 is considered. For the lon¬ 
gitudinal resonance, model I provides the correct qual¬ 
itative behavior and reproduces the local maximum ac¬ 
curate in position. In general, however, it does not pro¬ 
vide highly accurate quantitative agreement. Especially 


for larger anisotropies (wy/wz > 7), the model 1 over¬ 
estimates the coupling strengths. Again, this behavior 
is understandable since for the higher-order longitudinal 
resonances the anharmonicity becomes increasingly im¬ 
portant, which cannot be modeled accurately with wave- 
functions in the harmonic approximation. 

For the transversal resonances with A(o, 2 .o) 
A(2,o,o) c.m. excitation, the coupling strengths predicted 
by the model 1 agree quantitatively very well with the ab 
initio ones for LOyjiOz > 2. This agreement is remarkable, 
since no free parameters were used in the models. 


VI. 3D-2D TRANSITION 

In the following, the transition from a 3D to a quasi-2D 
confinement is considered. Again, as a realistic example 
two ^Li atoms in a sextic confinement of A=1000 nm and 
14 = 35 .OEi- are chosen. To obtain a pancake-shaped trap 
iOx and Wz are decreased, keeping the ratio = l-I 

constant. The lowest-order resonance with (0,2,0) c.m. 
excitation and the next to leading order resonances with 
n = (4,0,0) and n = (0,0,4) are considered. Again, in 
analogy to the ID case these resonances are the most 
pronounced ones, having in mind that the two lowest- 
order resonances with excitations in the weakly confined 
directions, i. e. with c.m. excitations (2, 0, 0) and (0, 0, 2), 
fade away towards large negative values of dy/a, loosing 
their resonant character. 



Figure 5. (Color online) Resonance positions for different c.m. 
excitations obtained by full ab initio calculations (dots) and 
the model (lines). The dashed lines are obtained by the per¬ 
turbative correction of the energies (model B), the solid lines 
indicate the numerically exact c.m. correction (model C) and 
the dashed-dotted lines correspond the harmonic approxima¬ 
tion (model A). 


As before, the behavior of the resonance positions in 
Fig. [5] is discussed first based on the ab initio results. 
A similar behavior as for the ID case is visible. By the 
same arguments that hold in the 3D to ID transition, 
the (0,2,0) resonance with excitation in the strongly 
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confined direction starts at negative values of dyja for 
small anisotropies and converges to an asymptotic value 
for strong anisotropies 0). The higher order 

(4,0,0) and (0,0,4) resonances start at positive values 
of dyja at small anisotropies and do not converge to an 
asymptotic value for uJx,‘^z 0. 

For the (0, 2, 0) resonance, the model A is again shifted 
to slightly higher resonance positions due to the ab¬ 
sence of a negative quartic term. While for intermedi¬ 
ate anisotropies the harmonic theory gives slightly bet¬ 
ter quantitative agreement to the ab initio calculations, 
the asymptotic value is quantitatively reproduced to high 
accuracy within the corrected model B, where a pertur¬ 
bative treatment of the energy corrections is sufficient. 

For the higher order (4,0,0) and (0,0,4) resonances, 
the perturbative treatment of the corrections gives an al¬ 
most perfect quantitative description of the resonance 
positions for mild anisotropies. However, for strong 
anisotropies it fails (shows a minimum in the resonance 
positions and then goes to positive values of dyja) and 
the exact treatment of the ID c.m. excitation within 
model C delivers the most accurate results. 



ment (that is why it cannot be resolved for 2), 

increases until it reaches a maximum and then falls off to 
zero as —>■ 0. Again, its behavior is a result of the 

counter-acting effect that on one hand decreasing the po¬ 
tential depth increases the anharmonicity and hence the 
coupling strength, but on the other hand —>■ 0 cor¬ 
responds to switching off the confinement leading to a 
vanishing confinement-induced coupling. 

A breakdown of model 2 is detected for very large 
anisotropies (wy/wz > 10) where the model 2 predicts 
negative coupling strengths for all resonances. The rea¬ 
son why the coupling integral in Eq. (22) can result in 
negative values is the negative quartic term in Eq. Q. 
However, while negative coupling strengths for them¬ 
selves are not a problem yet, an unphysical discontinuity 
is introduced when taking the absolute value. Hence, 
the sign change (or even vanishing value) of the coupling 
strength is unphysical. Since such a behavior is absent 
for the cigar-shaped regime, the used harmonic quasi-2D 
trap state wavefunction, Eq. (19) turns out to be inap¬ 
propriate here. 

Still, the model 2 reproduces correctly the decreasing 
coupling strength for the (0, 2, 0) resonance. Eor smaller 
anisotropies Wy/w^ < 5 it is even quantitatively accu¬ 
rate. For the (4,0,0) and (0, 0,4) resonances, the non¬ 
monotonic behavior is reproduced qualitatively. 

In general, for the positions as well as for the cou¬ 
pling strengths, the models in quasi ID show a better 
quantitative agreement than the corresponding models in 
quasi 2D, simply because for a single decreasing potential 
depth the anharmonicity effects are milder compared to 
the pancake-shaped potential and can be reproduced by 
the model that is based on the harmonic approximation 
more accurately. 


VII. SIMULTANEOUS VARIATION OF THE 
POTENTIAL DEPTH 


Figure 6. (Color online) Coupling strength for different c.m. 
excitations n = (n^,, riy, n^) (see legend) obtained by full 
ab initio calculations [dots (0,2,0), squares (0,0,4), triangles 
(4,0,0)] and the model 2 (solid lines). 

In Fig. the coupling strengths for the transition from 
a 3D to a quasi-2D confinement are shown. The ab initio 
results show a constant decrease for the (0,2,0) reso¬ 
nance. This is in analogy to the resonances with excita¬ 
tions in the strongly confined direction in the transition 
to a cigar-shaped potential shown in Fig. Again, the 
decrease can be explained by the decrease of the coupling 
potential W (r, R) cx 14 if 14 and 14 are reduced in the 
transition to a pancake-shaped confinement. 

The behavior of the coupling strengths of the (4, 0, 0) 
and (0, 0,4) resonances in a pancake-shaped confinement 
exhibits a similar behavior as the longitudinal (0, 0,4) 
resonance in a cigar-shaped confinement shown in Fig. 
The ab initio results demonstrate that the coupling 
strength is close to zero for an almost isotropic confine- 


It is an important question how the c.m.-rel. coupling 
can be optimized. As demonstrated above, the coupling 
at the lowest-order resonances, i. e. with c.m. excitations 
Hi = 2, = 0, and with excitations in the tightly 

confined direction have a peak coupling strength for an 
isotropic trap and then monotonically decrease with the 
anisotropy, i. e. with a decreasing potential depth in the 
weakly confined direction(s). 

Higher-order resonances in the weakly confined direc- 
tion(s) have a very small coupling for an isotropic con¬ 
finement, peak at mild anisotropies and then decrease to 
zero for increasing anisotropies. 

In general, the coupling strength can also be modified 
by a simultaneous variation of the potential depth, i. e. 
by a variation of all potential depths and not only se¬ 
lected ones leading to a different (quasi-lD or quasi-2D) 
trap geometry. The ab initio results in Table [H] demon¬ 
strate that the coupling increases with a decrease of the 
potential depth V. This behavior can be understood in- 
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Vy/Ey 

53.9 

36.0 

18.0 

9.0 

4.5 

initio [10 flLOy^ 

7.7 

9.3 

13.2 

19.0 

25.0 

Hniodel[10 ^hlUy] 

5.2 

5.9 

6.7 

6.3 

2.6 


Table II. Coupling strengths Wn for ujxjujy = 1, WyjuJz = 0.5 
of the (0, 2, 0) resonance for different values of the potential 
depth for ab initio calculation and model 1. 


tuitively. As the potential becomes deeper, the harmonic 
approximation becomes more accurate and it has a zero 
c.m.-rel. coupling. While the model follows this behavior 
for a deep potential, it looses the accuracy, if the poten¬ 
tial gets to shallow, and results in an unphysical decrease 
of the coupling. The reason for this failure of the model 
is that the harmonic wavefunctions become less accurate 
for a decreasing potential depth. 

In the case of an optical lattice, the decrease of the 
potential depth has yet another important consequence. 
A tunnel coupling between neighboring wells enhances 
drastically the c.m.-rel. coupling. In fact, ab initio calcu¬ 
lations of a double and a quadruple-well potential have 
shown that in this case even high-order resonances show a 
considerable coupling. Moreover, the results of the calcu¬ 
lation explain some of recently measured loss resonances 
in a shallow 3D optical lattice in an ultracold gas of Cs 
atoms [32] ■ 


Vy/Er 

53.9 

36.0 

18.0 

9.0 

4.5 

dy /Q-ab initio 

0.9 

0.91 

0.92 

0.9 

0.81 

dy j Q-rnodel 

0.82 

0.8 

0.77 

0.72 

0.66 


Table III. Resonance position for /ujy = 1, ujyjuj^ = 2 
of the (0, 2, 0) resonance for different values of the potential 
depth for ab initio calculations and the model B with the 
perturbative energy correction. The harmonic approximation 
of the model A predicts a resonance position of dy/a = 0.886 
independent of V. 


With a decreasing potential depth the harmonic ap¬ 
proximation of the potential becomes inaccurate and the 
limitations of the introduced models become visible. In 
Fig. i a good agreement of the harmonic approxima¬ 
tion was visible for the positions, especially for small 
anisotropies. In the harmonic approximation, i. e. with¬ 
out the energy corrections in Eqs. (121 and (di, the po¬ 
sition of the resonances is independent of the potential 
depth. However, with decreasing potential depth, the 
anharmonic terms in the sextic potential start to have 
a significant influence already at energies of the lowest 
trap state and hence influence the resonance positions. 
In Table |III| the dependence of the position for a mild 
anisotropy is compared for the model B and ab initio cal¬ 
culations. While the ab initio results are almost constant 
for a deep potential, the resonance position decreases for 
a decreasing potential depth to small values. While the 
model A in the harmonic approximation is independent 


of the potential depth, the energy-corrected model B re¬ 
flects this decrease reasonably. 


VIII. WAVEFUNCTION ANALYSIS 


C.m.-rel. coupling resonances have been directly ob¬ 
served experimentally in a two-body system via coher¬ 
ent molecule formation [12] and indirectly in a many- 
body system in terms of particle loss and heating urn- 
The latter is a consequence of the molecule formation at 
the resonance. Both measurements detected resonances 
where the bound state was excited in a strongly confined 
direction, i. e. in quasi-lD the (2,0,0) and (0,2,0) reso¬ 
nances jlldl], and in quasi-2D the (0,2,0) resonance |7] 
was detected. Due to the anisotropy of the confinements, 
the resonance position occurred for positive values of the 
s-wave scattering length. It was demonstrated above, 
see, e. g. Fig.|^ that also resonances at negative values of 
the s-wave scattering length occur: for lowest-order reso¬ 
nances if the anisotropy is kept small, or for higher-order 
resonances for a stronger anisotropy of the confinement. 

In the following, ab initio wavefunctions are analyzed 
for the system of two ^Li atoms confined to a sextic 
potential with parameters A=1000nm, = 1-1, 

oJxl^y = 10, Vy = 35.9 Ey. The corresponding energy 
spectrum is shown in Fig.j^ Considered are the densities 
of the bound and trap states involved in the transversally 
excited (2, 0, 0) resonance and the longitudinally excited 
(0, 0,4) resonance. The positions at which the wavefunc¬ 
tions are investigated are chosen such that the overlaps of 
the involved trap and the bound states are still small in 
order to compare the characteristics of the states. This 
is dy/a = 1.42 for the (2,0,0) and dy/a = —1.2 for the 
(0, 0,4) resonance, respectively. 

In Fig. cuts through the trap-state densities are 
shown. Since both states have the same diabatic state, 
i. e. |'i/'i‘J^(o,o.o)); they have the same global nodal struc¬ 
ture, i. e. two regions of large probability to the find the 
particles separated from each other, away from the diag¬ 
onal zi = Z 2 ■ This can also be seen considering the mean 
radial density 

nOO 

f= drrp(r). (24) 

Jo 


which is a measure for the mean distance of the particles. 


In Fq. (24), 


p{r) = ■ 


dFRdDr|^(r,R)|' 


(25) 


is the radial pair density where 4'(r,R) denotes the full 
six-dimensional wavefunction of the system, dVR is the 
c.m. volume element and dDr is the angular volume ele¬ 
ment of the rel. motion. For the trap states of Fig. I^the 
radial pair density is shown in Fig. The large probabil¬ 
ity for the particles to be off-diagonal in Fig. [^are clearly 
reflected. This can be quantified by the mean radial dis¬ 
tance which is ri = 3.95 dy = 1.25dz at dy/a = 1.42 and 
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Figure 7. (Color online) Cuts along the elongated z-direction 
= X 2 = 2/1 = 2/2 =0)1^) through the full six¬ 
dimensional ab initio wavefunction. 

Upper part: The ground trap state diabatically described by 
IV’i'^“(o,o,o)) at dy/a = 1.42. 

Lower part: The ground trap state diabatically described by 
IV'i'l'(0,0,0)) at dy/a = -1.2. 

(For both plots an identical color code has been used). 


¥/ = 4.56 dy = 1.44^2 at dy/a = —1.2. Hence, the mean 
distance of the trap state is on the order of the longitudi¬ 
nal trap length dz = 291 nm which reflects the elongated 
shape of the confinement. 

In the region of interaction a strong and small-scale 
nodal structure is visible close to the diagonal zi = 
where both particles are close to each other. The nodal 
structure is also visible in the inset of the radial pair- 
density plot which shows p{r) for small r. In this region, 
the Born-Oppenheimer interaction potential possesses a 
deep minimum (compared to the energy scale of the trap¬ 
ping potential) which supports many bound states lead¬ 
ing to the rich nodal structure. 

In Fig. 1^ cuts through the bound-state densities are 
shown. In both states, the particles only occupy regions 
where they are very close to each other, i. e. close to the 
diagonal Zi = Z 2 - The bound state at dy/a = 1.42 (upper 
panel in Fig.j^ has no c.m. excitation in the z direction. 



Figure 8. (Color online) Radial pair density of the hrst trap 
state for different values of the s-wave scattering length. 

Hence, no zeros (nodes) are visible in the density (wave- 
function) in scales of the trap length. Of course, the many 
small-scale oscillations in the bound-state regime stem¬ 
ming from the deep Born-Oppenheimer interaction po¬ 
tential are still present. The bound state at dy/a = —1.2 
shows four large-scale nodes along the z direction which 
is due to the (0,0,4) c.m. excitation of this bound state. 

At the transversally excited resonance, at dy/a = 1.42, 
see Fig. the atoms in the bound state which can be 
approximated by \4’{b) ‘J’( 2 ,o,o)) have a mean distance of 
rb = 0.29 dy, i. e. it is small compared to the confinement 
length in the tight y direction. This demonstrates the 
strong binding of the atoms. Away from the (0, 0,4) res¬ 
onance, at dy/a = —1.2, the atoms in the bound state 
have a mean distance of TT = 1-01 dy = 92.9 nm, i. e. it 
is on the order of the trap length in the tightly confined 
direction. 

An interesting question is whether the bound state at 
the (0,0,4) resonance at dy/a = —1.2 has enough bound 
character to trigger molecule formation and subsequent 
losses in a many-body system. 

In |12] a molecule formation was observed experimen¬ 
tally at the resonance where the atoms in the bound state 
had a mean radial distance of nf = 140 nm. This is even 
larger than the value oir/ = 96.8 nm at the (0,0,4) res¬ 
onance. Hence, a molecule formation with subsequent 
processes is also expected at this resonance. 

IX. CONCLUSION 

Experiments [a [12 have demonstrated that inelastic 
confinement induced-resonances can influence the stabil¬ 
ity of an ultracold quantum gas and can be adopted to 
create molecules fully coherently. The resonances de¬ 
tected so far were measured in a strongly anisotropic 
confinement at large positive values of the s-wave scat¬ 
tering length. In fact, it has been demonstrated that 
they were all of a special kind, namely the ones excited 
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Figure 9. (Color online) Cuts along the elongated z-direction 
= X 2 = 2/1 = 2/2 =0)1^) through the full six¬ 
dimensional ab initio wavefunction. 

Upper part: The c.m. excited bound state diabatically de¬ 
scribed by |'!/’ 6 'F( 2 ,o,o)) at dy/a = 1.42. 

Lower part: The c.m. excited bound state diabatically de¬ 
scribed by |V’ 6 ‘F(o,o, 4 )) at dy/a = —1.2. 

(For both plots an identical color code has been used.) 


in the strongly confined direction. In the present work it 
is demonstrated that also resonances in the weakly con¬ 
fined direction occur. Models are introduced to describe 
the resonance position and the coupling strength of, in 
principle, all c.m.-rel. coupling resonances of a system of 
two ultracold atoms. 

A study of the most pronounced resonances is per¬ 
formed for a variation of the external confinement. The 
lowest-order resonance in the strongly confined direc- 
tion(s) and the next to leading order resonances in the 
weakly confined direction(s) are discussed in the tran¬ 
sition of an almost isotropic 3D confinement to a quasi- 
ID (cigar-shaped) and a quasi-2D (pancake-shaped) con¬ 
finement. While the position and the coupling strength 
of the resonance(s) excited in the strongly confined di¬ 
rection converge monotonically to a constant (non-zero) 
value for an increasing anisotropy, the position of the res- 
onance(s) excited in the weakly confined direction fade 



r [a.u.] 


Figure 10. (Color online) Radial density of the c.m. excited 
bound states. 


away to negative infinity for an increasing anisotropy and 
the coupling strength approaches zero. These resonances 
show a maximum in the coupling strength for intermedi¬ 
ate anisotropies. 

The models are discussed in comparison to ab initio 
calculations. In the transition to a cigar-shaped con¬ 
finement geometry, the resonance positions are described 
accurately by the model C. The coupling strength is de¬ 
scribed quantitatively correct by model 1 for the reso¬ 
nances in the strongly confined direction and qualita¬ 
tively correct for the resonance in the weakly confined 
direction. 

In the transition to a pancake-shaped geometry, the 
resonance position of the resonance with an excitation in 
the strongly conhned direction is described quantitatively 
accurate. For the resonances in the weakly confined di¬ 
rection, the accuracy for different anisotropies depends 
on the used approximations. The coupling strengths are 
reproduced qualitatively by the model 2 except for large 
anisotropies. For the latter the model 2 faces limitations 
and results in negative coupling strengths. 

A variation of the potential depth shows that the c.m.- 
rel. coupling can be increased by decreasing the potential 
depth. While the position of the resonances in a very 
shallow sextic potential is still described accurately by 
the model B, the values for the coupling strengths loose 
their accuracy, since the harmonic approximation of the 
wavefunctions the model is based on looses its validity. 
Therefore, while the discussed models provide a helpful 
guide, the full ab initio calculation remains indispensable 
for highly precise quantitative predictions or for describ¬ 
ing properly some trap geometries. 

The analysis of the wavefunctions involved in the res¬ 
onances in a cigar-shaped potential demonstrates that 
molecule formation and subsequent losses are also ex¬ 
pected for the resonance excited in the weakly confined 
direction. In this case they occur for large negative val¬ 
ues of the s-wave scattering length. We hope that this 
type of resonance will soon be verified experimentally. 
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The study of inelastic CIR has demonstrated that one 
of the most fundamental and routinely adopted approx¬ 
imations in ultracold atomic quantum gases — the har¬ 
monic approximation — has to be abandoned in order to 
describe particle loss, heating, and molecule formation 
in a variety of experiments. In fact, the inelastic CIR 
can be tuned not only by a variation of the scattering 
length but alternatively by a modification of the confine¬ 
ment geometry. Hence, it might deliver a novel tool for 
ultracold-atom experiments to alter the interaction be¬ 
havior by a variation of the external confinement in the 
vicinity of an inelastic CIR. This can be valuable in cases 
where the standard technique of using magnetic Feshbach 


resonances may not be available, such as for earth-alkali 
atoms. Similar to a magnetic Feshbach resonance, at the 
inelastic CIR the association of molecules can be tuned 
fully coherently. 
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X. APPENDIX 


A. 2-Channel model 

To obtain the resonance position and coupling strength 
from ab initio calculations, a two-channel model is fitted 
to the eigenenergy spectrum in the vicinity of the investi¬ 
gated resonance. There are two diabatic states, the trap 
state |t) and the bound state \b) with diabatic energies Et 
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and Eb. Introducing a coupling W between these states, 
the Hamiltonian matrix 


(Et 

VVI^ Eb) 


(26) 


is obtained. A diagonalization of this matrix by a linear 
transformation = U~^E[U, where U consists of the 
eigenvectors of the diagonal matrix H^, leads to the en¬ 
ergies El and E 2 of the adiabatic states which are known 
from the ab initio calculations. Assuming that the dia- 
batic states are linear in the vicinity of the avoided cross¬ 
ing, i. e. Et = ax+ b and Eb = cx + d (where x = ^), the 
coefficients a, b, c, d and the coupling W can be obtained 
by a minimization \\U~^HU — Hd||- The position of the 
resonance is then easily obtained from the crossing point 
of Et{x) and Eb{x). 


in where M = 2m. For simplicity, only a single 

spatial direction is considered in the following and the 
subscript j is omitted. The first-order energy correction 
is determined by 

/ OO 

dR\'tlj{R)\^V''’^'>{R). (28) 

-00 

An exact expression for the integral of a triple product 
of Hermite polynomials and a Gaussian is known to be 

m 


dxe Hfc(x)H„(x)H„(x) = 


" y/irklnh 


(s — k)!(s — n)!(s — m)\ 
(29) 


where 2s = n + k + m must be even. To make use of 
this formula the and i?® terms in need to be 
expressed in Hermite polynomials. For example, 


B. Perturbation Theory 

For a correction of the resonance positions in the model 
B, the energies of the c.m. trap states in the sextic po¬ 
tential are treated within first-order perturbation theory. 
Since the c.m. Hamiltonian separates, it is sufficient to 
evaluate the ID Hamiltonians. The unperturbed system 
is the ID harmonic oscillator for which the wavefunctions 


are given in Eq. (15). The anharmonic terms of the sextic 
potential 


i?4 = lH4(A) + ^H2(i?) + 5Ho(i?). (30) 

Inserting the expressions for R'^ and i?® into the integral 
in Eq. (28), splitting the integrals, and evaluating each 


with formula Eq. (29) yields 


^(a) ^ _ 


1 


36 {2n^ + 2n + l)VH^uj^ 


1152 H2 

— (4 -I- 6 -I- 8 n -I- 3) 


(31) 




1 Hujj 


Rl 


720 H/ 




(27) 


are treated as a perturbation. Here the potential is writ¬ 
ten in dimensionless units of energies in Hujj and lengths 


To determine the 3D perturbative energies = 

y z of fflo sextic potential the anhar- 

fa) 

monic energy corrections E^ j of the three spatial direc¬ 
tions and the corresponding harmonic oscillator energies 
E^t'^ = Hujj {n+ need to be added up. 












